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Abstract. A mechanism deriving new well-posed evolutionary equations from 
given ones is inspected. It turns out that there is one particular spatial operator 
from which many of the standard evolutionary problems of mathematical physics 
can be generated by this abstract mechanism using suitable projections. The 
complexity of the dynamics of the phenomena considered can be described in 
terms of suitable material laws. The idea is illustrated with a number of concrete 
examples. 
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Introduction 

In [9j [12] it has been shown that the standard (autonomous and linear) initial boundary value 
problems of mathematical physics share a simple common form, if considered as first order 
systems. Indeed, it is found that they are of the form 

d V + AU = F, 

where in the usual cases A is skew-selfadjoint, and U and V are linked by a so-called material 
law 

V = M (d^ 1 ) 

where M is a bounded operator- valued function, analytic in a ball Be (r, r) of radius r e M>o 
centered at r. A4 (6q ) is then well-defined in terms of an operator- valued function calculus 
associated with do as a normal operator in the weighted L 2 -type space H gj o(M.,H), g > j-, 
with inner product 



(U,V)^ f (U(t)\V(t)) H exp(-2gt) dt. 

Jr 



We will not need to recall the solution theory of such equations, (which we like to refer to as 
"evolutionary" as the term "evolution equations" appears to be reserved for a rather special 
case in this wider class), since the purpose of this paper is not on well-posedness issues but on a 
remarkable even more specific structural similarity between various equations of mathematical 
physics. We shall indeed see, that most of the standard initial boundary value problems of 
mathematical physics can be derived from a single spatial differential operator of the form 

A =\y o 

with a suitable domain to make A skew-selfadjoint in 

h-.= [®L2(n))e(0^(n)). 

\fceN / \feeN / 

Here Q is a non-empty open subset of a Riemannian C^i-manifold M with V = d<£) denoting 
the co- variant derivative (Riemannian connection). The spaces 

are the completion of (real- or) complex-valued Lipschitz continuous covariant fc-tensor fields 
having compact support in Q with considered in the norm | • \ k induced by the inner product 



o,vo 



JM 



where V denotes the volume element associated with the Riemannian metric tensor field g 
given by the Riemannian structure of the manifold M. Here (ip\ili) k abbreviates the function 



where 

(*>*)" (*\*)k,(TM p )* 

is the (real) inner product of covariant fc-tensors on the tangent space TM p atpeM and 7TT 
denotes complex conjugation. Covariant 0-tensors are simply real numbers and so 0-tensors 
fields are real-valued functions on M and so we let 

($>*)^<$l*>0,(TM p )* :=$ *- 

Since for k e N>o covariant fc-tensors on the tangent space TM p are elements in the (real) 
tensor product space (X) fe (TM p ) the inner product is induced by 

<$ <8> • • • <8> $fc-l|*0 ® • • • ® ®k-l>k,(TM p )* = <^o|*0>(TA/ p )* • • • (®k-l\Vk-l) (T M p )* ■ 

In the sense of this inner product V* = — div is the formal adjoint of the co- variant derivative 

V. 

The complexity of the various physical phenomena has no influence on the choice of A but 
is reflected in different material lawajj. The process of extraction of particular operators 
from the "mother" operator (BqM (Bq ) + A) is surprisingly simply and amount to projecting 
this operator (class) down to smaller subspaces. The transparency and simplicity of this 
construction is rather striking and shows the interconnectedness of various different physical 
phenomena if inspected from a mathematical point of view. This connection becomes obscured 
if a second order (or higher order) model is chosen as a starting point. It appears that a largely 
misguided pre-occupation with the occurrence of the Laplaciaro in equations of mathematical 
physics has lead to a dominance of second-order equations and systems in various models. As 
it turns out, however, the first order approach leads to a unified and transparent access to a 
large class (if not all) of typical linear model equations. 

In the applications we shall for sake of simplicity focus on the Cartesian or periodic case 
(M = R n - k x T fc , n = 1,2,3, k = 0,1,2,3, k ^ n, with T being the flat torus obtained 
from the unit interval [—1/2, l/2[ by "gluing" the end points together (implying periodicity 
boundary conditions on ] — 1/2, l/2[ in the last k components), which is perfectly sufficient to 
understand the reduction mechanism and its applicability. 



1 Mother Operator, Relatives and Descendants 

1.1 A Construction Mechanism 

Definition 1.1. Let C : D (C) Q Hq — > H\ linear, closed and densely defined and B : Hq —* X 
a continuous linear mapping, X, Hq, Hi Hilbert spaces. We say B is compatible with C if 



^^This is the opposite point of view to the "conservation law" perspective. 

2 This is surely fostered by the comforting regularity properties of elliptic (and parabolic) differential operators. 
In our perspective these have their place when qualitative properties are of prominent interest, but, as it 
turns out, have limited importance and are often distractive for fundamental well-posedness issues. 
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• CB* is densely denned (in X). 

Theorem 1.2. Let C : D (C) E Hq — » H\ linear, closed and densely defined and B : Hq — > X 
a continuous linear mapping, X, Hq, Hi Hilbert spaces. Moreover, let B be compatible with 
C. Then 

(CB*)* = BC*. 

Proof. It is 

CB* c (BC*)* . 

Let now ue D ((BC*)*) then 

(v\(BC*)*u) = (BC*v\u) 
= (C*v\B*u). 



We read off that 

and 

Consequently, we have 

and so 



B*ueD(C**) = D(C) 

CB*u= (BC*)*u. 

CB* = (BC*)* 



(CB*)* = (BC*)** = BC* 



□ 



Definition 1.3. Let C : D (C) Q Hq — ► H\ linear, closed and densely defined, X, Hq, H\ 
Hilbert spaces. Moreover, let Bq : Hq — > X be compatible with C and B\ : H\ — > Y be 
compatible with C* . Then we call B\CBq the (Bq, !?i)-relative (or simply a relative) of C. If 
not both of the mappings Bq, B± are bijections, then we call B\CBq the (.Bo>I?i)-clescendant 
(or simply a descendant) of C (and C the mother operator of B\CBq). 

Of particular interest are compatible operators resulting from orthogonal projectors. We 
introduce the following schemes of notation: 

Definition 1.4. Let V be a closed subspace of a Hilbert space H. Then we denote the 
orthogonal projector onto V by Py. 



Definition 1.5. Let X®Y a direct sum of Hilbert spaces X, Y. Then we denote the canonical 
projections 

X@Y — X , X@Y ^Y 
x®y i-» x , x®y i-* y 

by ttx, t^y, respectively. 

These notations are employed in the following elementary observation, which we record without 
giving the elementary proof. 



1.2 Evolutionary Operators and their Relatives 
Proposition 1.6. Let V be a closed subspace of a Hilbert space H . Then 

Py = TTyTTy. 

and 

(1 - Py) = ■Ky ± TTy± = Py±. 

Moreover, 7iy is the canonical isometric embedding ofV in H and 

TTyTTy and Try±TTyj_ 

are the identities on V and V^, respectively. 

The subspace V ©{0} of V ®V is commonly identified with V. We shall, however, prefer to 
distinguish Py and iry, since the latter allows a proper formulation of reducing an operator 
equation to a subspace by constructing appropriate descendants. 



1.2 Evolutionary Operators and their Relatives 

The concepts introduced in the previous section extends naturally to evolutionary problems 
as described in the introduction. 

To be specific we consider a particular class of evolutionary problems of the form: 



{d M{d 1 )+A)U = F (1) 

where 

A -\c 

and C : D (C) E Hq — > H\ is a closed densely defined linear operator so that A is skew- 
selfadjoint in the Hilbert space H ■■= Hq H\. We assume for the material law that there is 
a Co £ M>o with 

** (x R<0 (™o) U\d M (do" 1 ) U) > co <x R<0 (mo) U\U) (2) 

for all sufficiently large g e M>o and all U £ D (do) E H e (R, H). This assumption warrants 
solvability and causality of the solution operator, for this variant of the solution theory compare 

In order to ensure ^ we may and will assume that 

M (d, 1 ) =M + tfM® (d^) 

with Mo selfadjoint and ^ m [H]Mo^ Mo[hv ^[{o^Ma ^ (M m ( S o~ 1 )) f [{o)]M uniformly strictly 
positive definite for all sufficiently large g e M>o- If 

M™ (d, 1 ) = Mi + M® (do 1 ) (3) 
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with the operator norm in L (H gfl (M, [{0}] Mq) , H gfi (M, [{0}] Mq)) 

K[{0}]Mo ^ (Mi) 7T[*o}]Mo) ^{OPo ^ e ^ <2) (^o" 1 ) 7r f{0}]A4 < X ( 4 ) 

then a standard perturbation argument shows that (J2J) is maintained. Since solution theory is 
not the topic of this paper we will not dwell on these issues in the following. We merely note 
that the construction of relatives and descendants maintains this solvability condition. 

We note first that as a by-product of the above we have the following. 

Corollary 1.7. Let C : D(C) ^ Hq — * H\ linear, closed and densely defined, and let Bq : 
Hq — > X be compatible with C and B\ : H\ —* Y be compatible with C* . Then Bq © B\ is 
compatible with 

' -C* 

Moreover, we have that the relative 



A:- 



(S © B X ) A (So © B{f = [b^ b * BqC q * BI ) (5) 

of A is skew-self adjoint in X © Y . 
Proof. Note that 

*«*-(?*)-(?!)(U)-(U)(?! 

Then the equation follows from the previous theorem by noting that Bq compatible with C im- 
plies So©l is compatible with A and similarly l©Si is compatible with (Bq © 1) A (Bq © 1)* = 

I since B\ compatible with C* implies B\ compatible with BqC* = (CSq)*. 
CB Q U J 

Indeed, since 

B C*Bl c B~q~C*BI 

the density of the domain follows. □ 

-C 

C 

and (BqM (d^ ) — A) are relatives via I I as a unitary mapping. 

Applying earlier observations to evolutionary operators yields the following result. 

Theorem 1.9. Let C : D (C) c Hq — » Hi linear, closed and densely defined, and let Bq : 
Hq — > X be compatible with C and B\ : H\ — ► Y be compatible with C* . Moreover, let 

1 C 



Remark 1.8. The structure of A as I _, I implies that the operators (BqM (Bq ) + A) 



Then (Bq © B\) A (Bq ® B\) is skew-self adjoint and 



Bq (Bq © Si) M V (Bq © Si)* + S © B X A (Bq © Si 



is an evolutionary relative of (BqM (d Q ) + A) . 



1.3 Reducing Equations to Subspaces 

Unitary equivalence of A would be a typical example illustrating the previous theorem. More 
interesting, however, are descendants being produced by projections onto proper subspaces. 



Remark 1.10. In contrast to the last result general relatives of (poM (8q ) + A) need not 
maintain its clear formal structure. Indeed, a few elementary row and and column operations 
can produce almost impenetrably confusing "model equations". In many instances it turns out 
to be the main task to reconstruct (poM (8q ) + A) from a quite different looking (often only 
a formal) relative. 

1.3 Reducing Equations to Subspaces 
1.3.1 The General Construction 

Let now V be a closed subspace such that 

P V M (c^ 1 ) = M [d l ) P V . 
Then, applying iry to equation (|13p we obtain similarly as before 

8q (nyM (df) 1 ) 7ly) KyU + 7Ty A (llylTyU + TTyxTTyXU) = TTyF. (6) 

This is now an equation for ixyu in V with little chance of being well-posed (due to the free 
floating part ir v± iry±u). To enforce well-posedness, we could assume that iry and Tr v are 
compatible with A and instead (assuming tt v u = 0) consider 

(do (nyM (Oq 1 ) tt v ) + (vrv/Avr^)) ir v u = tt v F. 



In general, however, isyAiXy = iryAir v will fail to be skew-selfadjoint although 

inherits its positive definiteness property from A4 (<9q~ ) : 
There is a Co 6 M>o with 

&* (x R<0 (mo) U\do7ryM (Bq 1 ) tt v u) n >c (x (m ) tt v u\tt v u) 

\ / £)]0,0 \ / g,0,0 



co{X R<0 {m )u\u 

' e,0,0 



(7) 



for all sufficiently large q e M>o and all u e D (do) E H e (R, V). 

But if V = VoSVq, i.e. iry = iry ©7Tvi, where Try is compatible with C and iry 1 is compatible 
with C*, we see that we have an evolutionary descendant. We summarize this observation in 
our next theorem. 

Theorem 1.11. Let Vo E Ho, V\ E H\ he closed subspaces such that iry is compatible 
with C and iry 1 is compatible with C* . Then with V ■= Vo © V± we have that the operator 

(do (jryM (do 1 ) n v ) + {^yA-Ky) ) is the (^,71^) —descendant of (poM (d^ 1 ) + A) . 
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1.3.2 A Particular Case: Removing Null Spaces 

To see the above construction at work let us consider the standard issue of reducing A to the 
ortho-complement of its kernel for the operator 

(d M{d^)+A). 

We assume for sake of definiteness that 

M {d, 1 ) =M + d Q l M x (do" 1 ) . 

We have 



*-m M ( 5 ° *) ^^i) + ^mm^ 



A[H] 

and 

(do (nmiAM (Bq 1 ) 7rf {0}] 



are two corresponding relatives. Note that here ^ A v H i A is already closed and vrrjojiA 7 ! = 0. 
Moreover, 



VT-T777T = TTT^TrTFT ttf VT 7 



A[H] ~ "C*[#i] ^ "C[H ] 

and 

^[{O}]^ = 7r [{0}]C ©7T[{0}]C*. 

Clearly, vrrfQijc and ^ r<*\H 1 are com P a tible with C and correspondingly 7rr/o}]c* and 7 r„r„ -, are 
compatible with C*. So, if at least one of the null spaces [{0}] C or [{0}] C* is non-trivial, which 
is the only interesting case, we have that (do (^~X[htM {^o ) ^Trm ) + ^ATflT^TTrm" ) an< ^ 

( <?0 (vr[{o}]A-^ (^o~ ) ''"[Wi.a) ) are m( ieed descendants of (poM (<9g ) + A) . How can these 
descendants help in solving a problem for the mother operator? 

To simplify calculations we first confirm that we may assume that Mq can be replaced by 
^MolH^MolH] =Pm [H]- Indeed, 

H = Mo [H] © [{0}] Mo 

M = ^M [H]Mo-KM o[H] ®O m]Mo 

Mo = ^M [H]M Tr* Mo[H] ®TT[{0}]M ^[{0}]Mo 

we obtain 



and 



With 



M l Mo\IM Q 1 = \/M l \jMo {1& M h] © [{0}]A4o ) (l Mo[H] © O m]Mo ) \JMo\/M 1 
= (l«b[fl] © °[{0}]M ) = ^om^oM- 



1.3 Reducing Equations to Subspaces 



By writing again A for J M~ Q 1 A^M 1 and M (d x ) for J M 1 M (d x ) J M 1 ([I]) is re- 
covered but now with a material law operator of the form 

Since we have 

H = A[H]®[{0}]A 

we obtain the decomposition 

[ d (^^mM (d Q X ) vr^) + tt^Att^-,) vr- 



TU + 



'A[H] JVI \ v J "A\W\J ^ "A[H]~"A[H]J "A[H]' 
+ d ° ( 7T A[H] M (^) *[{0y\A) 7T [{0}]A U = ^A[H] F 
and 

[So (n m]A M (d^ 1 ) vrf {0}]A )) i [{0 pu + (do (n [{0}]A M fc 1 ) ^y)) k^u = n m]A F. 
Note that 



7P 



Att^-, 



A[H]~ n A[H] 



is now skew-selfadjoint in A [H] and strict positive defmiteness of the real parts of the operators 

^° ( f i[m ^ (^o" 1 ) ^"Trm") anc ^ ^° ( ^[MJA-M (^o" 1 ) ^UoWA ) * s P re served. It should also be clear 
that here it makes no sense to assume vrrjoin^u = since the system would be over-determined 
since in general we cannot assume that -P[{o}]A commutes with M (d^ ). Instead, solving the 
latter equation for 7Tr/o}]A w yields 

K[{o}]AU = [Tr m]A M (Bq 1 ) irf {0}]A ) d Q ~ 1 7r[ { o}] A F+ (8) 

- (n m]A M (Oq 1 ) Trf {0}]A ) (v m]A M (d^ 1 ) n^j ^^sf. (9) 

Inserting this into the first equation yields 

(d (kj^M (d^) vr^) + n^An^) vr^ 



t u+ 



d o ( 7r A[F] M ( 5 ° ^ 7r [*o}]A) (^[{0}]A-^ (d X ) tf { o }]A ) (vr [{0}]A 7W (d p x 



^Ajti}) "A[fl]' 



= (doM (do 1 ) + tt^Att^) IT^U (10) 

Now, if 



M (d^) ■■= TTAmM fa 1 ) 7T* 



^Aim M Vo 1 ) ^{oha) (^[{o}]aM {d^ 1 ) rf mA ) * (kuokaM fa 1 ) Tr*^ 



fill 



satisfies the required strictly positive defmiteness for material laws in the usual sense in the 
subspace A [H] , we have solvability in A [H] . This observation for general A may be used to 
restrict the problem to A [H] , which may allow for example to utilize the compactness of the 



2 Some Applications 



restricted resolvent of A in A [H] (if this indeed ho lds) whereas the original A may have too 
large a kernel to have a compact resolvent, see e.g. |16j for an application. 

( - \ t - N- 1 / - \ ( 12 ) 

- {^A[m doM ( cl o ) "ko)]a) ^[{0}]a3 X (5 : ) n? {0}]X J [7T m]A d M (8 J ) -Tr^^j 

Assuming Q, (HJ), it suffices to inspect (|lip for .M (5^ ) replaced by .Mo + do Ml due to the 
smallness assumption ([!]). Indeed, for the regular case where Mq is strictly positive, according 
to the above, it suffices to consider M.$ = 1 for which 



^{ u \ 8 ^Am 7r Am u / e , 



for ue H eA [R,A[H]\. 

Substituting the solution of this standard evolutionary problem into ([S]) we obtain the other 
part of u. 



2 Some Applications 

2.1 A Particular Mother Operator 

Now we consider specifically 

{d M t^ 1 ) + A) U = F (13) 

where 

'-(ST) 

with a suitable domain making A skew-self adjoint in the Hilbert spacqfl 

ff=(®i*(n))©(e4(n)). 

\fceN / \fceN / 

Here V and V* are formal adjoints on the linear subspace 

(0(7 ljfc (o))e(©c' 1)fc (o)) 

VfceN / VfceN / 



3 We have chosen here to add up tensor spaces of all ranks, although in applications only k = 0, 1, 2, 3 appear 
to be relevant. Restricting to the "physically relevant" subspace 

V = L 2 (Q) ® L\ (fi) © L% (Q) © h\ (fi) 

may therefore be considered as a first application of the mechanism to generate descendants of the operator 
in JHTU). 
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2.2 Isolated Physical Phenomena 

of H, where C\^ (fl) denotes the space of Ci-smooth co-variant tensor fields of rank k with 
compact support on a Riemannian C^i-manifold M with metric tensor g. The differential 
operator V is the so-called co-variant derivative and its skew-adjoint —V* is frequently in- 
troduced as the tensorial divergence^] div. For sake of definiteness we shall only consider the 
choice^ 

( - ( V" ' 



A:=\ 

V v o y 

where V denotes the closure of V applied to elements of ©fceN ^lk (^) as an operator in 
©fceN-^fc (O) . For the material law we impose the usual constraint ([2|). 

It may be surprising that the majority of initial boundary value problems from classical math- 
ematical physics can be produced precisely from (|13ll4p by choosing suitable projections for 
constructing descendants. This is the main application of the above considerations. In order 
to make matters more easily digestible we constrain the illustration of our observations to the 
simple flat case, i.e. M is W 1 ^ x T fc , n = 1,2,3, k = 0,1,2,3, k < n, where T is the flat 
Torus. 

It will turn out that the physical interpretation has remarkably little relevance for our struc- 
tural observation. In fact, problems can be very different in physical interpretation, sharing 
the same formal structure makes the solution theory coincide. 

2.2 Isolated Physical Phenomena 

2.2.1 Acoustic Equation, Heat Conduction and the Relativistic Schrodinger Equations 

If we choose 

v = l§ (n) e £? (n) , n^m 3 , 

for our construction of descendants via projectors then we obtain the classical system governing 
acoustic waves or, merely depending on the choice of material law, the heat equation. In the 
Cartesian case we have by identifying 0— tensors with functions and 1— tensors with vector 
fields the classical first order system 



with 



M(tf) = ('° K ) + *Z 



pv\ , a -i f oo 

a 



as a simple material law operator. If p, k are strictly positive definite, continuous, selfadjoint 
operators then the material law can be taken to describe acoustic wave propagation. If a 
is also strictly positive definite, continuous and selfadjoint operators, we interpret a as a 
damping term. Alternatively this could then be considered as describing heat propagation 



4 Correspondingly, we could use grad as a notation for the covariant derivative V to give the original evolu- 
tionary equation the suggestive look of the acoustic system (see below) . 
5 This choice may be referred to as the Dirichlet boundary condition case. 
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2 Some Applications 



with Cattaneo modification. Keeping all these constraints except for assuming k = 0, we get 
the classical heat propagation. The second row describes in this interpretation (for g = 0) the 
so-called Fourier law of heat conduction. In both cases, the materials are indeed such that 
the material law commutes with complex conjugation. This allows to interpret the equation 
in real- valued terms. 

Alternatively, we may view L? (Q) = L 2 , (fi,C), k £ N, as a Hilbert space over the field K. (by 
restricting the underlying scalar field). Then we have 



TZ:Ll(n,C) 



u 



Ll(n,w)@L 2 k (n,: 

9teu' 

3m u 



as an M-unitary mapping. For example multiplication by the complex unit is then unitarily 
equivalent to 

1 



1I1K- 



With this observation we get that the Schrodinger operator do + i Ajj assumes the real form 

*+(aVo D )' < 16 > 

where, to be specific about boundary conditions, we have chosen the Dirichlet-Laplacian Ajj. 
Clearly, due to its second order type this operator is not covered in our approach. There is, 
however, a variant known as the relativistic Schrodinger operator in which —Ap is simply 

, which turns out to be essentially unitarily equivalent to the 



replaces by ^/—Ao 



grad 



acoustics problem (|15|) for a even simpler material law. 



Indeed, removing the null space of A 



div 
grad 



by reducing further to the subspace 



L (fi) © grad [L 2 (fi)], which is the range of A, we obtain another descendant of (|13ll4p . 
which is indeed a relative of the relativistic Schrodinger operator. According to the polar 
decomposition theorem there is a unitary mapping U such that 



grad = U grad 



and 



Consequently, 



div 



grad 



U*. 




M(d { 



d M {d { 
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2.2 Isolated Physical Phenomena 



Since for the relativistic Schrodinger operator M. (d ) = 1, we get 



d + 



div 
grad 



1 
U 



do + 




grad 



grad 




1 

U* 



There is another rather common version to translate the operator of the wave equation into 
a first order in time system, which , however, is nothing but another relative of the acoustics 
operator. Utilizing the naive analogy to the ordinary differential equations case one translates 



BqU - A D u = f 



into the system 

d + 



£ 





0A o -e 
1 



d u 
—u 



d + 



A D 

1 



d u 

—u 



Here we choose — e £ M^o m the resolvent set of —Ap. The reasoning goes like this: 

0A„-£ 
1 

is skew-selfadjoint considered in L 2 (J7)©iJi (V - ^D + e) where Hi U/—Ad + e) is D (V~ A~d) 
equipped with the inner product 



(u,v) »->■ (J-A D u\^-A D v\ +e(u\v) , 



i.e. for e = 1 the graph inner product of \/—A]j. We use that 

V-Ad + e : Hi (y/-A D + e) - L 2 (fi) 



0A D 
1 



i I e\\ , fOA D -£ 



is unitary. Transferring this considerations to a more general situation we note 
d M (V) - 

= a °(W> +a »~ i too)) + (r D o~ 



and 







An-e\ fl 



y/-A D +eJ \l 








^-A D + e~ 



_ -V-A D + e 
V-A D + e 

1 n( . ° - 

/ \ grad + iy^ 



grad 



VM/i o 
u+ 
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2 Some Applications 



Note V-Ad + e 



implies 



grad 



We 



using polar decomposition for 



grad 



We 



grad 



We = U+ 



u+ 



grad 
grad 



IV e 



grad 



grad 



iv e 



grad 



iVe 



ive 



we 



grad 



We 



and 



t/| = t/+. 



With the polar decomposition of grad = U 



grad 



U :=U®h 



Note 
where 



U 



grad 



we define 

[{0}]gr°ad- 

[L 2 (0)] - grad [L 2 (fi)] , 



grad 



[L 2 (0)] = ([{0}] 



grad 



[{0}]grad)" 
div [L 2 (0)] . 



With this we get 




We 



div +iy/EU 

grad + iy^^ 

U* 




div \ 
grad ) +W£ {U 

Thus we obtain for the transformed equation a new material law: 

1 \ . . /__1x /l 



A4 (d^) 



UU-y/-A D + e 



-m a- 1 ; 







-l, 



y/-A D + e l U+U* 



-l, 



-i / eV-Ao + e L U+U* + U/eU 



iy/eU 







Given the complexity of the arguments needed it is still a surprisingly common mechanism 
used (at least in the case e = and simple material laws) to turn partial differential equations 
of wave equation type into first-order-in-time systems. None the less, in the above terminology 
we encounter here mere relatives of the acoustic equations, which in turn is a descendant of 
(fT3TT4l) . 
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2.2 Isolated Physical Phenomena 

2.2.2 Elastic Waves 

If we choose 

V = L\ (0) © sym \L\ (SI)] 

with Q open and non-empty in R 3 then we obtain the classical system governing the prop- 
agation of waves in elastic or, depending on the choice of material law, viscoelastic waves. 
Here sym is the mapping sym : L 2 ($7) — > L 2 , ($7) induced by the symmetrization operation for 
co- variant tensors of rank 2 

Thus, in Cartesian coordinates the elasticity operator is 

where 

Gradu = - (diVj + 5 i i; i ) iJ=lj2i 3 

and its negative adjoint 

Div T = div T = ( J] c^- J 

Vi =1 / i=l,2,3 

for suitable displacement velocities t> = (vj) •_ x 2 3 an d symmetric stress tensors (Tij)^ . =1 2 3 . 
A discussion of various possible material laws of interest can be found in [9] , compare [I] . 

2.2.3 Electro-Magnetic Waves 

If we choose 

V = L\ (0) © asym [L| (ft)] 

with ft open and non-empty in IR 3 then we obtain the classical system governing the propa- 
gation of electro-magnetic waves. Here asym is the mapping asym : L\ (ft) —* L\ (ft) induced 
by the anti-symmetrization operation for co- variant tensors of rank 2 

(%, y) •-* g ( T ( x > y) ~ T (?/' x )) 

To convince ourselves that this leads to Maxwell's equations we calculate V • W for anti- 
symmetric tensors W £ Coo, 2 (ft) m Cartesian coordinates, ft being a non-empty open subset 
ofM 3 . Let© 

W = oJie 2 © e 3 + w 2 e 3 © e 1 + uj 3 e l © e 2 - wie 3 © e 2 - a^e 1 © e 3 - w 3 e 2 © e 1 

then 



TU+i = ^fe+2e ®e 



fc+i 



W = asym (T ks e k ® e s ) = I (T fea e fe <g> e s - T fcs e a ® e fe ) = T fe3 e fe a e 3 = 2 ^ fc<s mod 3 T ka e k a e s 
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2 Some Applications 



V-W = d k W ks e s 

= d]Wi2e 2 + d!W 13 e 3 + d 2 W 2 3e 3 + d 2 W 2 ie 1 + d 3 W 32 e 2 + d 3 W 31 e l 
= d\u 3 e 2 — d\oj 2 e 3 + d 2 ui\e 3 — d 2 u) 3 e — d 3 uj\e 2 + d 3 u 2 e l 
= (d\uj 3 - d 3 uji) e 2 + (d 2 0Ji - d\0J 2 ) e 3 + (d 3 oj 2 - d 2 u 3 ) e 1 



curl | <jj 2 | . 
w 3 . 



Correspondingly, 

V = r/ie 1 + i] 2 e 2 + r/ 3 e 3 



and so 



W = dxrue 1 (g) e 1 + di^e 1 <g) e 2 + dir^e 1 (g) e 3 + <? 2 r/ie 2 (g> e x + 

+ d 2 r] 2 e 2 ® e 2 + <? 2 r/ 3 e 2 (g) e 3 + <9 3 r/ie 3 (g) e 1 + <9 3 r/ 2 e 3 (g) e 2 + d 3 r/ 3 e 3 (g) e 3 . 

We see that 

asym (VV) = 
= d x n 2 - (e 1 ® e 2 - e 2 (g) e 1 ) + <9ir/ 3 - (e 1 (g) e 3 - e 3 (g) e 1 ) + 5 2 ??i- (e 2 (g) e 1 - e 1 (g) e 2 ) 

+ d 2 j] 3 - (e 2 (g) e 3 - e 3 (g) e 2 ) + <9 3 7/i- (e 3 (g) e 1 - e 1 (g) e 3 ) + d 3 r] 2 - (e 3 (g) e 2 - e 2 (g) e 3 ) 

= (<3i?72 - d 2 r?i) 2 ( el ® e 2 - e 2 (g) e 1 ) + (<9 2 r/ 3 - d 3 n 2 ) - (e 2 (g) e 3 - e 3 (g) e 2 ) + 

+ (5 3 r?i - 5x77s) g ( e3 ® el - el ® e 3 ) 

= (dim - d 2 rji) e 1 a e 2 + (d 2 n 3 - d 3 n 2 ) e 2 a e 3 + (<9 3 r7i - dir/ 3 ) e 3 a e 1 

= d A 7] 



and so 



asymV —■ cIa 



on differentiable 1-form fields. In Cartesian coordinates Maxwell's equations assume the fa- 
miliar form 

curl 



*m & ') ■ , cuil „ 



where 



curl = curl 



and containment of £" in D I curl J encodes and generalizes the electric boundary condition, i.e. 
vanishing of the tangential components of E, for the electric field E to the arbitrary boundary 
off). 
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2.2 Isolated Physical Phenomena 

2.2.4 Reducing Dimensions 

Another instance of the reduction procedure under discussion is the reduction of the dimension. 
Indeed, the subspace V to which the reduction of (|13[ [bi)) could be discussed, might be 
described as the null space [{0}] B of a suitable operator B (or a family of such). An example 
in case would be to take the last partial derivative d n in a periodic direction and consider 
V k := [{0}] d n c L 2 k (Q) , k e N. In other words, we assume that Sl = l!oxTc R™- 1 x T. Of 
course there are many other options of this kind for example looking at radially independent 
solutions in polar coordinates etc.. In order to keep matters simple we shall here only consider 
the case M = R n_1 x T. In this case, e.g. 7iy is compatible with V and ttv 1 is compatible with 

I V I . Consequently, the reduction process can be utilized to reduce the spatial dimension by 
one. 

Similarly, in the case = Q,q x T n_1 E R x T" -1 = M (implying periodicity in the last (n — 1) 
coordinates), we may reduce the problem to one spatial dimension by taking the projection 
onto V ■= fl?=2 [{0}] dj- Resulting in a (1 + 1) -dimensional evolutionary descendant of (|13|14p 
on the open subset £1 of the flat tube manifold M. We may write this descendant in Cartesian 
coordinates simply as 

d M(d ^ + ^y a?) 

8 A : . _i ir.ji.^ ._ .1 r2/n \^ T2 



where now I ° I is skew-selfadjoint on the space L (P-o) © L (Oo) • 

\di ) 

If XIq = K we may go one step further, we can decompose L 2 (R) into orthogonal subspaces 

L 2 (R) = L 2 ' even (R) © L 2 ' odd (R) , 



with 



L 2,even ^ = {f e L 2 (R) \f (x) = f (-x) for a.e. x 6 R} , 
L 2,odd ^ = {feL 2 (R) |/ (x) = -/ (-x) for a.e. x e R} 



Since TT L 2,eveiim\ and ir L2 ,eveii( R \ are compatible with d\ = d\ we obtain that 



d x 
d\ 



2,even m\ m r 2,odd 



is the f7r L 2,even (R ),vr odd J -descendant of JT7J on L 2 ' even (R) © L 



with 



V L2.odd (R) y y u 7r L 2,odd ( 

as a new material law operator. If Ai (d^ ) is block diagonal 

M(d-^-( Mo0 ^ ° 



17 



2 Some Applications 

then the two rows can be combined into one 



d i^T L 2,even m M o {d X ) < 2 ,even (K) + 7r L2: odd (R) - M n ( 5 o 



')':,odd (l 



on L 2 (R) = L 2 ' even (R)©L 2 ' odd (R) . This is the so-called transport equation in the (1 + 1) -di- 
mensional case, which thus also is shown to be a descendant of (|13|14p (for fi = M = lx T n ). 

Remark 2.1. (The Transport Equation in R n ) Returning to R n and assuming that by a suitable 
choice of coordinates the transport operator 8q + a ■ d assumes the unitarily equivalent form 

do + dx 

with d\ on a cylinder^] fi := R x fio c R x R n_1 . Here also d\ = d\. Of course, we could have 
more complicated material laws: 

5 X(5 - 1 )+5 1 . (18) 

The cross-section Q,q of the cylinder serves here merely as a parameter range, since no differ- 
entiations in these directions are involved. Allowing for additional parameter dependence in 
(|13|14p would make JTHJ a descendant of (|13|14p . 

2.3 Interacting Descendants 

The various descendants of (|13|14p can interact in many ways to create new models of more 
complex phenomena. We shall first discuss a particular interaction based on an alternating dif- 
ferential forms framework (alternating covariant tensors) . Then we shall turn to the discussion 
of coupled descendants, where the coupling only occurs via the material law operator. 

2.3.1 The Extended Maxwell System and the Dirac Equation 

The Extended Maxwell Operator Assuming a relatively simple material law of the form 

M {d Q l ) = Mo 

with A^o continuous, selfadjoint and strictly positive definite, Maxwell's equations can be 
reformulated as 



\diA / 

Here d\ a is the exterior derivative applied to covariant 1-tensors (with Dirichlet type boundary 
condition). By including alternating tensor fields of all odd orders in the first block component 
and of all even orders in the second block component we arrive at 



'o -i- v m | v / | v m 

dl )3 A 



^Similarly we may consider transport on a period slab, i.e. SJ = Tx!]o cTxRasa flat Riemannian manifold. 
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2.3 Interacting Descendants 

Here di^Ais the exterior derivative applied to the direct sum of alternating tensors of order 1 
and 3 (with Dirichlet type boundary condition). Note that d a lj = on 3- forms in M 3 . The 
material law operator Mq is here of course assumed to be continuous, selfadjoint and strictly 
positive definite on the larger space. Adding 



-d ,2A 




VMi 



as another descendant^ of 



-V* 
V 



we obtain a unitarily equivalent variant of the type of operator discussed in [8] as the extended 
Maxwell operator 




-4,2 a 




y/Mo- 



Here do 2 A is the exterior derivative applied to the direct sum of alternating tensors of order 
and 2 (with Dirichlet type boundary condition). In Cartesian coordinates this is 



Mo" 1 



/0 \ 
- curl 


\ curl J 



M r 



VM 



/ div\ 

grad 

div 

Vgrad / 



Vm 



(19) 



which is a convenient reformulation of Maxwell's equations for regularity and numerical pur- 
poses, compare |15| . That the spatial part of this extended system is still skew-selfadjoint and 
that it can be reduced to the original Maxwell system is due to the fact that 



Mr 



/o \ 

- curl 

\ curl J 



M^\ y/Mr 



/ div\ 

grad 

div 

\grad / 



VMj 



are commuting selfadjoint operators, which are indeed annihilating each other. The possibility 
of reconstructing the original Maxwell system assumes a particular form_| of the right-hand 
side, see |8j for details. By adding a material law term M. [p^ ) we can allow for more 
complicated material behavior: 



^Note that 



—do, 2 a 

* 
do, 2 a ) 



-1 

1 





do, 2 A 



do, 2 A 





1 
-1 



For these special data it can be shown that components of order and 3 are actually zero. If general right- 
hand sides are considered then these components will be non-zero producing what is called "scalar waves" 
contributions. 
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2 Some Applications 



d + M {d^) + 



Mo 1 



/O \ 
- curl 
OOOO 

\ curl 0/ 



M^ + VMo 



/ div\ 

grad 

div 

\grad 0/ 



Vm, 



Remark 2.2. Projecting this further down by eliminating the third row and column leads to 
a slightly smaller descendant of the extended Maxwell system 



d + M (d 



-v 







I -curl 
curl 



div 


grad 



For "ellipticizing" Maxwell's equations, e.g. for numerical purposes, this modification is per- 
fectly sufficient, [T7] , 



The Dirac Operator The Dirac operator Qo{do,d) is usually given as the (4 x 4)— partial 
differential expression with the block matrix form (for mass equal to 1) 



Qo(d ,d) 



(do + i C(8)' 
\C(d) d -i 



Herg^|C(d) := 



<9 3 dt-i di 
<5i + id 2 — d 3 



Sfe =1 n fc 4, where 



°-i?iV' 



-i\ ._ +1 

+i o ' ns - -1 



J Note that 



d 3 d\—i d 2 

8i +182 -d 3 



is an operator quaternion, since it has the form 



A -B* 
B A* 



where A : D (A) c H —* H, B : D (B) c H — > H are closed densely defined linear operators, such that A 
has a non-empty resolvent set g {A) and A, B* axe commuting, i.e. 

{X-Ay^cBiX-Ay 1 

for A e g(A). If A,B are complex numbers (as multipliers) this block operator matrix yields a standard 
representation of the classical quaternions. 
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2.3 Interacting Descendants 



are known as Pauli matrices. Applying the unitary transformation given by the block matrix 
j^f+i +1 
V2 \ i -1 



to Qo (^0 , d) we obtain 



Qi(d ,d) :- 



do i-iC(d) 
i + iC(d) d 

\-i +1 



i -1 

-i +1 
i -1 



-i do + 1 + C(d) 
-i C(d) + d - i 

do + i C(d) 
C(d) do - i 



-id + l-C(d) 
-i C(d) -d + i 




The latter may be a preferable form since Qi(do,d) has the typical Hamiltonian form of 
reversibly evolutionary expressions of mathematical physics 

do -W* 
W d 



where W 



iC(d) 



ids + i ic?i + di 

id\ — di — ids + i 

On first glance the Dirac operator does not seem to fit into the framework we are discussing 
here, since it does not appear to be constructed from descendants of (|13|14|) . A closer inspec- 
tion, however, shows that the Dirac operator is actually unitarily equivalent to, i.e. in the 
above sense a relative of , the extended Maxwell operator (with a variant of the material law) . 

To see this connection we separate real and imaginary parts, which yields that W corresponds 
to 

\ /0 -1 \ 

10 

0-1 



/ 



1 



V 





+ d 3 
-d 2 

di 



■1 - d 3 d 2 
di 

-di 

-di 1 - 5 3 



-di 

di 

-l + d 3 





\0 1 / 



ds 
-d 2 



-d 3 



-di 

-d 2 



di 



-d 3 



-d,\ 
di 
d 3 




Noting that 

/ 1 

0-10 
0-1 
-10 



0\( 

d 3 
-d 2 



o / 



\dl 



-d 3 d 2 

di d 2 

-d t d 3 

-d 2 -d 3 / 



di\ / 1 \ 

10 
10 

\o i oj 



/ o ft d 2 d 3 \ 

-d\ d 3 -di 

-di -d 3 d\ 

\-d 3 di -d l / 

div 
grad curl 



and 



/ 1 0\/0-lfl 0\ 



0-10 
0-1 
-10 / 



10 

0-1 

\o 1 J 



/0 1\ 
10 
10 

\o i o/ 



/ 1 \ 

0-10 

0-1 

V-i 0/ 

/0 1\ 

10 

0-100 
\l o o/ 



/-l 0\ 
1 
0-10 

V o i o o/ 



21 



2 Some Applications 



we obtain the unitary equivalence 



/ /0 1\ 

10 

10 

\0 1 0/ 

/0 0\ 





V \0 0/ 

/ 



/0 0\ 



\o o o oj 

1 




-10 



1 



1 





w 



-w* 





//0100\ 

10 

1 

\1 0/ 

/O 0\ 





V \o o o oy 



div 
grad curl 



V 





In the free-space situation the Dirac operator 3q 
to the extended Maxwell operator 




W 



-w* 





is thus unitarily equivalent 



d + M± 



div \ 

grad — curl 
div 
\grad curl / 



where M\ 




is skew-selfadjoint, i.e. from 



the electrodynamics perspective we are in a chiral media case. 



Thus we have shown that the Dirac equation also fits seamlessly into our construction of 
descendants of ()13I14|) and their interaction. In particular, the Dirac operator is a relativd 11 ! 
of the extended Maxwell operator discussed above. 

Remark 2.3. It is a rather remarkable observation that the Dirac equation is so closely con- 
nected to the extended Maxwell system. It appears from this perspective that spinors are 
actually a redundant construction since the alternating forms setup for the extended Maxwell 



11 In the framework of quaternions a connection between a differently extended time-harmonic Maxwell op- 
erator and the time-harmonic Dirac operator has earlier been discovered by Kravchenko and Shapiro, [2], 
compare also [?]. 
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2.3 Interacting Descendants 

system is already quite sufficient to discuss Dirac equations. The interpretation of this ob- 
servation is not a mathematical issue but may well be a matter for theoretical physicists to 
contemplate. 

2.3.2 Coupled Systems 

Let us recall from |12] the systematic coupling mechanism between various different descen- 
dants. Without coupling the systems of interest can be combined simply by writing them 
together in diagonal block operator matrix form: 



fV \ 



\VnJ 



A 



fU \ 



\U n ) 



fk\ 



\fnj 



where 



f Aq ••• \ 

'■• : 

: '-. 

V ••• A n J 

inherits the skew-selfadjointness in H = ©fc=o n ^k from its skew-selfadjoint diagonal entries 
Ak : D (Ah) ^ Hk — * Hk, k = 0, . . . ,n. The combined material laws here take the simple 
block diagonal form 



V 



fV \ 



\V n J 



M in (Bq 1 ) u 



fM 00 (do l ) o 




\ fU \ 



M nn (Bq 1 ) J \u n J 



Coupling between these phenomena now can be modeled by expanding the material law to 
contain block off-diagonal entries 



/ Mqo (do 1 ) 



ex /a— 1> 



m cx (d 



M 0n (d 1 )\ fM 00 (d 1 ) 




\M n0 (d 1 ) M nn (d 1 )J V 

The full material law now is of the familiar form 

V = M (do 1 ) U 
with 

M (d l ) ■= M in a 1 ) + M cx (do 1 ) . 



\ 



M nn (d Q - 1 )) 
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2 Some Applications 

This coupling mechanism now allows to model thermo-elasticity, thermo-piezo-electro-magnetic 
and so on. A number of examples for coupled systems have been discussed elsewhere, see 
|13| El [5j [TTJ [IUJ [6] . The "philosophy" of this coupling mechanism is that coupling occurs only 
via the material law. 

In the following we shall illustrate the abstract coupling mechanism with a particular concrete 
example, which will at the same time serve to exemplify the construction of descendants of 
coupled systems. 

Starting point of our example collection is the classical system of thermo-elasticity, which will 
also allow us to re-iterate the point made previously in connection with the Dirac operator and 
the extended Maxwell operator that indeed systems with very different physical interpretations 
may share the same solution theory with differences being incorporated merely in possibly 
different material laws (or their interpretation). 



Thermo-Elasticity and Biot's Model for Porous Media The classical system of (1 + 3) —di- 
mensional thermo-elasticitjo can be described by 



{d M (do" 1 ) + A) 



fv\ 



\Tj 



with 



A 



I -div \ 

-grad 

-Div 

\ -Grad / 



The classical material law is of the form 
with 



Mr 



( gi + r*c- 1 r 








r * c -l\ 




















82 





V c- l v 








c- 1 J 



M 



i •= 









.-1 

v 
















0/ 



where Qi,k,Q2,C continuous selfadjoint and strictly positive definite. This system formally 
coincides with Biot's porous media model, merely the meaning of the quantities involved, i.e. 
the units, have changed, see e.g. [5]- 



Due to an inconvenient choice of unknowns the original classical system of thermo-elasticity has an un- 
bounded coupling term, see |3j. The form given here avoids this drawback. More specifically the difference 
hinges on the use of stress instead of strain as unknown tensor field. 
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2.3 Interacting Descendants 



As in all these models we may allow for more complex material laws as long as ([2]) is maintained: 



/ 



d M (d- 



V 





t 





—div 








o-V 




-grad 












-Div 




{ 








-Grad 






fv\ 



\tJ 



F. 



(20) 



Reissner-Mindlin Plate Assuming li:=(l xTcR 2 xT=:M (instead of M = M 3 ) we can 
reduce (|20|) by one spatial dimension to a (l+2)-dimensional evolutionary problem following 
the strategy in Section 12.2.41 Indeed, the resulting evolutionary equation looks the same, but 
now it has to be interpreted in h\ (fi ) © h\ (f2 ) © h\ (fl ) © sym \L\ (fio)] with U Q B 2 . 
With 

ff\ 




and 



9 

w 



M (d l ) =M + d l Mi 



with 



f 61 








^ 




(d 








°\ 






K 






02 






, Mi-.= 








1 


-1 








\o 








c- 1 ) 




\o 








0/ 



Mr 



where Qi,k,Q2,C continuous selfadjoint and strictly positive definite (with physically different 
meaning, i.e. different units!) we obtain the Reissner-Mindlin plate model. Coupling occurs 
here via A4\. 

Note that by reducing this to a second order system (by substituting the equations from rows 
2 and 4 into the remaining two equations) we obtain the perhaps more familiar form of the 
Reissner-Mindlin model (with homogeneous Dirichlet boundary condition) 



QiSqTJ — div k x ( grad?] + s J + ddorj 
62SqS — Div CGrads + k~ 1 I gradr; + s 



/■ 
(]■ 



(21) 



where rj ■■= d n, s ■■= d s. 

For the damping coefficient d = we have that the system is conservative, since 



/ 



Mr 



/0 0\ 



0-10 
10 

\\o o o o/ 



/ -div ° \\ 

grad 

-Div 

-Grad 



V 



/ 



Mr 



is skew-selfadjoint and thus generates a unitary group leading to norm conservation for pure 
initial value problems. 
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Remark 2.4. 

1. (A note on the Kirchhoff-Love plate) 
Letting in 



( 0l 








^ 




/d 








°\ 





K 
















-1 











Qi 





T" 





1 








\o 








c- 1 ) 




\o 








0/ 



.4 



k = and £2 = (in consequence destroying well-posedness for associated initial bound- 
ary value problems) yields 



do 



o 



Vo 


















\ 




c- 1 J 



(d 








0/ 



A 



Eliminating the second and third unknowns and equations yields 
d 



QX 

o c- 1 



d 




div Div 

Grad grad 



This is the Kirchhoff-Love plate model, which by a suitable choice of boundary condition 
is again accessible to the abstract solution theory of evolutionary equations, see e.g. |12| . 



Love plate we could also formally obtain the real Schrddinger operator do - 



2. Following the "logic" of the transition from the Reissner-Mindlin plate to the Kirchoff- 

-A c 

.Ad 
see Section I2.2.1| from the first order system 



/1000\ /0 0\ 
OeOO 00-10 

OOeO + 0100 

\0 1/ \o o o o/ 



V 




'-Ad 


-V-Ad 
















V-a d 







\ 



/ 



by similarly letting e 
tions. 



V^Ad 
and eliminating the second and third components and equa- 



The Timoshenko Beam Assuming S\ 



fii x T 



x T =: M (instead of Q c 



for the Reissner-Mindlin plate, following the arguments in Section 12.2.41 we can reduce this 
model further to the (1 + 1) -dimensional case, which leads to the Timoshenko beam model. 
In Cartesian coordinates this is now 



/ 



d M (d t 



-i> 



V 



( -di 
-5i 
-3i 
V -bx / 



\\f v \ 

c 



where the material law has the same shape as before 



\Tj 





9 

\0j 
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with 



M {do 1 ) =M + d Q l M x 



References 



Mo 








K 






\ 












02 






c- 1 ) 



, Mi 



Id 








°\ 








-1 








1 








\o 








0/ 



where qi,k,Q2,C continuous selfadjoint and strictly positive definite. Of course physical 

meaning (the units) have changed once again. Here also the second order system may be more 

familiar 

gid^rj - diK' 1 (diTJ+a) + dd rj = f, 

(22) 



Q2%s — d\Cd{s + K (dirj+s 



9, 



where rj •■= d Q ?], s ■= d Q s, compare ([21D . 

Remark 2.5. (Euler-Bernoulli Beam) Repeating the questionable "construction" of the Kirchhoff- 
Love plate model for the Timoshenko beam, leads to the Euler-Bernoulli beam model. 

With this we conclude our tour through various examples underscoring the deep connected- 
ness of seemingly very different mathematical models. We have seen, how various particular 
dynamic linear model equations can be extracted from the mother operator (J13I14P assuming 
different material laws. The solution theory itself rests simply on strict positive definiteness. 
One may well wonder if the simplicity and transparency of these structural observations could 
not give rise to a "grand unified" numerical scheme. 
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